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Abstract 

The new parameterization of form factors developed for 4ir channels of the r 
lepton decay and based on Novosibirsk data on e + e~ — ► Att has been coded in 
a form suitable for the TAUOLA Monte Carlo package. Comparison with results 
from TAUOLA using another parameterization, i.e. the CLEO version of 1998 is also 
included. 



To be submitted to Comput. Phys. Commun. 



t Work supported in part by the Polish State Committee grants KBN 5P03B09320, 5P03B10721 and 
NATO Grant PST.CLG.977751, European Commission 5-th Framework contract HPRN-CT-2000-00149. 



1 Introduction 



Hadronic decays of the r lepton offer a unique laboratory for studying hadronic interac- 
tions at low energies (below the r lepton mass M). Proper modeling of such processes will 
facilitate comparison of experiment and theoretical models and may thus offer a crucial 
hint toward better understanding of low energy phenomenology of strong interactions. In 
this work we consider the r decay into four pions and neutrino. Several models of this 
decay are at present in use [1-6]. While Refs. [1,2] make an attempt to construct the 
Lagrangian of the decay from general theoretical principles, other models are purely phe- 
nomenological using either the experimental information from r decays like in Refs. [3,4] 
or from e + e~ annihilation like in Refs. [5,6]. In the following let us concentrate on the 
latter approach. 

Production of four pions is one of the dominant processes of e + e~ annihilation into 
hadrons in the energy range from 1.05 to 2.5 GeV. The hypothesis of the conserved 
vector current (CVC) relates to each other the cross section of this process and r — > 4ttv 
decay [7]. Therefore, all realistic models which describe the first process should also 
properly describe the other one. For a recent review of theoretical predictions for various 
decay modes of the r based on CVC see e.g. [8]. 

For more than five years, two new detectors CMD-2 [9] and SND [10] have been 
studying low energy e + e~ annihilation at the e + e" collider VEPP-2M at Novosibirsk. 
Their results can be used to provide a new parameterization of the form factors used in 
the Monte Carlo generators of the r lepton. In the present paper we define the form 
factors which can be used as one of the possibilities in TAUOLA [11-13], the r decay 
library. The TAUOLA package is organized in such a way that the phase space generation 
and calculation of the electroweak part of the matrix element are separated from the 
part of the code calculating the hadronic current provided by the particular model. This 
is convenient not only for constructing the Monte Carlo program, but also facilitates 
comparisons between different models. The technical description of the solution which 
allows an easy replacement of different hadronic currents is described in [14,15]. 

Our paper is organized as follows: in the next section we recall basic principles of how 
e + e~ data can be used in defining hadronic currents for the r lepton decay. In Section 
3 we describe those parts of the TAUOLA Monte Carlo algorithm for four pion generation 
and general formalism for the semileptonic/semihadronic decays which are needed for 
the definition of hadronic currents. In Section 4 we describe the new current for the 
47r decay of the r while in Section 5 details of the functions and constants used in this 
work are presented. Section 6 is devoted to technical tests of our code. Section 7 shows 
the comparison of numerical results from running the new version of TAUOLA with the 
relatively old one based on the CLEO parameterization [16-19] which is still widely in 
use despite the recent significant progress achieved by CLEO [3,4]. Chapter 8 briefly 
summarizes the main results of the paper. 

In Appendices we describe some technical aspects of our model. In particular, we 
tabulate the functions G(Q 2 ) used for the definition of the hadronic current. 
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2 Relation between r decays and e + e annihilation 
cross sections 

Via the hypothesis of conserved vector current one can relate the charged vector current 
coupling to the 4tt system to the electromagnetic (neutral vector) current measured by 
a(e + e~ — > 7 — > 4tt). There are two possible final states in e+e~ annihilation 1 

+ - * ~o -- + + 

e e ^7 — > p — > 7T 7T TT TT , 

e+e" - 7 * -> p° -> 7T + 7r-7r 7r°, (1) 

They are accessible by a different I3 component of the same 1 = 1 weak current describing 
r decay: 

r+ -> W+z/ T -> p+z/ T -> z/.tt+ttW, 

r+ -> W+z/ T -> p+P T -> z/ T 7T + 7r-7r + 7r . (2) 
The relations between processes (1) and (2) can be written as 

T(p+ -> tt+tt-tt+tt ) = ir(p -> tt+tt— tt+tt— ) + r(p° -> tt+tt-ttV), (3) 

T(p+ -> TT + TrW) = ^(p -> TT+TT-TT + TT-). (4) 

The prediction for T(r+ — > v T X + ) can thus be obtained from e+e - — > 4tt data and isospin 
invariance. Such a procedure was successfully applied in [5] where high statistics e+e~ 
data from the CMD-2 detector collected at center of mass energies (Q) from 1.05 to 1.38 
GeV were used [6]. Since e+e~ experiments are performed at fixed (at a time) Q, the 
integrated decay rates of p° to the 4-7T are well measured as a function of Q. Thus, the 
natural way of the r decay generation is to generate first the mass of the 4-7T system in 
accordance with the experimental distribution. Then, for the fixed Q 2 , the p° — > 4tt decay 
is generated. This is why in Ref. [5] the dT/dQ 2 distribution is generated independently of 
the differential distribution within the 4-7T system. The program is written and optimized 
to get maximum possible information from the experimental data. 

The approach of TAUOLA is somewhat different. A matrix element well isolated in a 
program module is separated into the hadronic and leptonic current. All physical as- 
sumptions on hadronic interactions are located in the hadronic current which features 
intermediate state resonances as well as other properties of the hard process. The phase 
space density generation and in particular the appropriate jacobians for mapping random 
numbers to phase space coordinates, which are defined independently from the particular 
process, are also calculated in the separate part of the code. Such an approach provides 
flexibility in studying particular choices of hadronic currents. In addition, the Q 2 distri- 
bution is not an input, but originates from the (partial) Monte Carlo integration of the 
matrix element over phase space. 

The aim of this paper is to use the model of [5] in the approach of TAUOLA which is 
more natural for comparison with other theoretical calculations/models. 

1 The e+e" -> 7r 7r°7r 7r° channel is forbidden by isospin and charge conjugation invariance. 
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3 General formalism for semileptonic decays 



The matrix element used in TAUOLA for the semileptonic decay t(P,s) — > u T (N)X is 
written in the form: 

M = ^=u(N)Y(v + a l5 )u(P)J^ (5) 

where J M =< X\V P — A^O > denotes the matrix element of the V — A current, relevant 
for the specific final state X. In general, the current J M depends on the momenta of all 
hadrons. N and P denote the four-momenta of the v T and r respectively. The squared 
matrix element for the decay of r with mass M and spin s reads: 

\M\ 2 = G 2 ^±^(c + ^), 



w = p"(n /t - 7f)0 ny, 

1 



H » = - p , pu )(K - 7«.n„) (6) 



with 



n M = 2[(j* • n) j, + (j • ao j; - {J* ■ J)i\g, 

IT 5 " = 2 Im e^ pa J* v J p N a , 

(7„ a = 1 in the Standard Model). If a more general coupling i> + 075 for the r current and 
v T mass m,/ 7^ are expected to be used, one has to add the following terms to u and H^. 

Co = 2t^^m u M{J* ■ J), 
v 2 + a 2 

2 2 

H» = -2 V -^- 2 m v Im e^.P u J p P a . (8) 
v 2 + a 2 

To obtain the polarimeter vector h in the r rest frame, it is sufficient to calculate the 
space components of h p = (H p + H p )/{uj + Cj) and set ho = 0. The differential partial 
width for the channel under consideration reads: 

2 2 

dT x = G 2 ^^rfLips(P; q h N)(cu + cu+ (H, + Hjs"). (9) 

The phase space distribution for the final state with four mesons plus neutrino is given 
by the following expression where a compact notation with q 5 = N and qf = mf is used, 

I 

dLips(P; q u q 2 , q 3 , g 4 , ?s) = T^Tf dC ? dC & dC & x 
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AP dQ * dQs Ql dQ2 Ql 



where 



Q 2 = (9i + Q2 + Q3 + g 4 ) 2 , Ql = (qi + Q2 + q 3 f, Ql = (?i + <? 2 ) 2 , 
Qmin = "ii + m 2 + m 3 + m 4 , Q ma;c = M - m 5 , 

<53,mm = mi + m 2 + m 3 , Qs.maz = Q - m 4 , 

<52,mm = ^1 + m 2 , <52,max = <?3 ~ ™3- (11) 

Here <if2 5 = dcosO^dn^ is the solid angle element of the momentum of v T in the rest 
frame of t(P), dVt^ = d cos 4 d7r 4 is the solid angle element of q± in the rest frame of 
Qi + <?2 + <?3 + ^3 = <icos#3<i7r3 is the solid angle element of <f 3 in the q% + q% + g 3 
rest frame, and finally, dVt 2 = dcos6 2 dir2 is the solid angle element of q 2 in the q% + g 2 
rest frame. 

These formula if used directly, are inefficient for a Monte Carlo algorithm if sharp peaks 
due to resonances in the intermediate states are present. We refer to the TAUOLA docu- 
mentation [11, 13] for details of the algorithm actually in use. For the present paper it is 
enough to note that those changes affect the program efficiency, but the actual density of 
the phase space remains intact. No approximations are introduced. 



4 Hadronic current for 4n system 

The model of Ref. [5] is based on the assumption that the ai(1260)7r and lutt intermediate 
states (which well describe the experiments on e + e~ — > 4n [6]), are dominant in the 
amplitudes r + — > v T f) + — > ^ r (4vr) + . In Ref. [5] it was shown that various two- and 
three-pion invariant mass distributions predicted by the model well describe experimental 
observations of CLEO [3] and ALEPH [20]. We include into consideration two most 
important channels of the a\ — > 3n decay (ai — > pn — > 3n and a\ — > an — > 3n) as well as 
the uj — > p7r — > 37r channel. Then for the process r + — > z/ T 7r + 7r°7r 7r° the current J M reads 

For the process r + — > z/ r 7r + 7r~7r + 7r°, where the a; meson also contributes, it is 

J M = J^-^pn + Jcn^ait + J^fm (13) 

where in the following we neglect the interference between the uj and ai currents. The 
z/ r 7r + 7r _ 7r + 7r° final states produced with the two currents (ai(1260)7r and ujtx) are effectively 
treated as distinct tau decay modes. 
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4.1 r + — > z/ r 7r + 7r°7r 7r° decay channel 

For the r + — > v t tt + (q^ir (q 2 )Tr° (q^Ti (q4) channel the current which includes possible 
Feynman diagrams 2 , can be written in the following way: 

(Q 2 )[ti (?2, 93, 9l, 54) + *i (?2, 94, 91, 93) + *i (93, 92, 91, 94) 

+ *i (93, 94, 9i, 92) + *i (94, 92, 9i, 9s) + *i (94, 9s, 9i, 92)], (14) 
Ja^air = G 7T+7T o 7T ^(Q 2 )[t^(q 2 , 91, 9s, 94) + 1%(qs, 91, 92, 94) + ^(94, 91, 93, 92) 

- *2 (9l, 92, 93, 94) - *2 (9l, 93, 92, 94) - ^2 (9l, 94, 93, 92)]- (15) 

Four- vectors t% and t£ have the following forms, where Q denotes Q = q\ + 92 + 93 + 94: 
*i (9i, 92, 93, 94 = n 77, T777 : r X 

DaAQ - qi) D p{<i3 + 94) 
{Q ■ (Q - qiM(Q - 91) • 93 - <g(Q - 91) • 9 4 ] 

+(Q" - rf)[(Q ■ 94)(9i • 93) - (Q • 93)(94 • 9i)]} (16) 
92, 93, 94) = = "/^'f^ S x 

DaAQ - qi) D Aq3 + 94) 

{q$Q • (Q - qi)(Q - 9i) 2 + (rf " Q")[(Q ' 9 2 )(Q - 9i) 2 ]}- (17) 

Here l/D ai (q), 1/D p (q) and 1/D a (q) are propagators of the a±, p and cr mesons, F ai (q) 
is the form factor and z is the dimensionless complex constant characterizing the relative 
fraction of the an intermediate state in the ai(1260) decay. G^+^o^o is some function 
depending on Q 2 which we find by fitting the 4n invariant mass distribution 3 . As a F ai 
form factor, we used the function from [5,6], F(q) = (l + m 2 ai / k 2 ) / {1 + q 2 / k 2 ) with A ~ 1 
GeV. 

The form of the propagators is very important for analyzing the data. We represent 
the function D{q) in the form used in [5,6] 



D{q) = q 2 -M 2 + iMT^L (18) 

g{M 2 ) 



where M and T are the mass and width of the corresponding particle, and the function 
g(s) describes the dependence of the width on virtuality. In the case of the p meson the 
function g p (s) reads: 

g p {s) = s~ l l 2 {s-4m 2 fl 2 , (19) 

while for the a meson it is: 

9_As) = (s - Arrf/s) 1 ' 2 , (20) 

2 Possible Feynman diagrams for r+ decays into 4ir via the aiir and the lot: intermediate states are 
shown in Appendix A. 

3 For details see chapter The G(Q 2 ) functions. 
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where m is the pion mass. 

The function g ai in the a\ propagator has the form: 

6 2 Pi ~ glP2 g2P3 ~ £3P2 ^y / SP2 

+ P2) D p (p 2 + p 3 ) D a (pi + p 3 ) 




\z\ 2 s 



3! 



Pi 



+ 



P2 



+ 



P3 



D a (p 2 +p 3 ) D a (p 1 +p 3 ) D a (p 1 +p 2 ) 



x 



x 



dpi dp 2 dp 3 <5 (4) (pi + p 2 + p 3 - 9) 



(21) 



2e 1 2e 2 2e 3 (27r) 5 

where q° = y/s, q = and Pi = (e,, p,) are the pion momenta in the rest frame of 
the p (the center of mass frame of the 4ir system). The first term corresponds to the 
ai — > 7r + 7r~7T decay while the second one to the ai — > 37r° decay. 



4.2 



-+ 



z/ r 7r + 7r 7T + 7r° decay channel 



For the r + — > z/ T 7r + (g 1 )7r (q , 2) 7r+ (?3)7T°(94) channel the current which includes the contri- 
bution from the ui meson intermediate state can be written in the following way: 



(Q 2 )[*l(9l> 92, 93, 94) + *i (?3, 92, 9l, 94) + *i (91, 93, 92, 94) 

+ *i (93, 9i, 92, 94) + *i (94, 93, 9i, 92) + *i (94, 9i, 93, 92)], 
W +XQ 2 ) [*2 (94, 9s, 9i, 92) + ^2 (94, 9i, 93, 92) 



Qi\ — XT7T 



TV 

" U)— >pTX 



- t 2 (qi, 94, 93, 92) - *2 (93, 94, 9i, 92)], 

G '^+7r- 7 r+ 7 ro( ( 5 2 )[*3(9l, 92, 93, 94) + t$(q 3 , 92, 9l, 94) ~ *3 (ft, 93, 92, 94) 

- t$(q 3 , 9i, 92, 94) - *s (9i, 94, 93, 92) - ^3 (93, 94, 9i, 52)]- 
Here four- vectors t%, t 2 are the same as in the previous case and t 3 reads: 



(22) 
(23) 
(24) 



t 3 (9l, 92, 93, 94) 



X 



A,(Q - 9i)^p(93 + 94) 

{92 [(Q ■ 9s) (91 • 94) - (Q ■ 94) (91 • 93)] - Q ■ 92(93 (91 • 94) - 94(91 • 9s)] 

+ (9i-92)[9 3 i (Q-94)- 94 l (Q-9 3 )]}, (25) 

where F u (q) and D w (q) are the form factor and propagator for the 00. Because of the small 
width of the u we set F w (q) = 1 and g u (s) = 1. G^+^-^+^o and G^+^-^+^o are functions 
of Q-. 
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Table 1: The masses and widths for the intermediate states. 



Intermediate 


Mass, GeV 


Width, GeV 


state 






p(770) 


0.7761 


0.1445 


ai(1260) 


1.23 


0.45 


u(782) 


0.782 


0.00841 


a 


0.8 


0.8 



5 Parameters used in the model 

The constants, widths and other parameters used in our numerical results were mainly 
taken from [21], for the p meson recent CMD-2 results were used [22]. Some of them are 
collected in Table 1, The m u = was assumed and PHOTOS [23,24] for QED radiative 
corrections was switched off. We used a r lepton mass of m T —\.lll GeV and the physical 
masses for the pions in the phase space, i.e. m n ± = 0.13957018 GeV, m n o = 0.1349766 
GeV. Note that because of CVC the masses of ^ and n° were set equal m^o = m^± = 
0.139570 in the hadronic current. The parameter A = 1.2 GeV for the F ai form factor 
was taken. The dimensionless complex constant z for t^qi, <?2, Q3, Q4) and for g ai (s) was 
set to (1.269, 0.591) which means that the ratio of contributions of the a to p intermediate 
states is 0.3 accordingly to Ref. [26]. 

5.1 The G{Q 2 ) functions 

New channels in TAU0LA are based on the same matrix elements for the p + decay as in [5] . 
The difference is in Q 2 generation. To complete the definition of the hadronic current, 
the appropriate choice of the three functions G(Q 2 ) was needed. We fit the functions 
G 7r + 7r o 7r o 7r o((5 2 ), G 7r + 7r - 7r + 7r o(<5 2 ), and G^ +/n .- 7V+7v0 (Q 2 ) in such a way, that TAU0LA reproduces 
predictions from the calculation [5] for the single variable dT/dQ 2 distribution. Samples 
of 5 000 000 events were used. The table of numerical values for the fitted functions is 
given in Appendix B. The functions are well determined by the e + e~ data from 0.9 to 
1.7 GeV, for lower energies the uncertainties are larger since the measured cross section 
is small, and we rely on assumptions made in Ref. [5]. 

Finally, we have normalized the functions G(Q 2 ) in such a way that our program 
reproduces r decay rates of Ref. [21], see results in Table 2. 
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Table 2: Numerical results for the integrated decay rates of our model 



Channel 


Decay rate, GeV 


T(r+ -> z7 r 7r°7r 7r u 7r+) 


2.4462* 10~ 14 


T(r + — > z/ r 7r + 7r + 7r~7r°) 


9.5130* 10- 14 



6 Numerical tests 

6.1 Technical test 

While developing the Monte Carlo algorithm, it is important to perform numerous tech- 
nical tests of the Monte Carlo code. Tests of the TAUOLA Monte Carlo library consisting 
of the comparison of numerical results obtained from the program and independent semi- 
analytical calculations are described and listed in Ref. [13]. At that time the agreement 
between the Monte Carlo results and some analytical calculations was pushed to the level 
of few permille. Recently, thanks to much faster computers, checks at the 0.1 % level 
were in some cases redone. In those tests different assumptions about matrix elements 
were used, e.g. masses of n^ 1 and it were set equal or even zero. The actual choice of 
the current was also modified. Parameters of the presamples were varied and results were 
checked to be independent of that. Satisfactory agreement was always found. This gives 
us confidence in the technical side of the algorithm. 

6.2 General test 

We have compared various possible invariant mass distributions constructed from the 
momenta of the r decay products in the Novosibirsk model coded in TAUOLA on one side 
to the code used in Ref. [5] on the other side. The samples used in the plots were of the 
size of the experimental data used as an input to the model, i.e. 30 000 events. In each 
case agreement was sufficiently good, for details see [25]. In this way we could convince 
ourselves that the adaptation was really correct. 

7 Comparison with other parameterizations 

Once we became confident in the correct implementation of the new Novosibirsk model, 
we have decided to compare its predictions to another one. For that purpose we have 
chosen the model tuned to the 1998 CLEO data because it is still quite often used. 

In this section we collect predictions for all possible invariant mass distributions for 
the decay channels: r — > z/ r 7r + 7r~7r + 7r° and r — > z/ r 7r + 7r 7r 7r . In Figs. 1 to 7 we compare 
predictions of our new Novosibirsk current with the older CLEO current [1, 17-19], both 
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Figure 1: The z/ r 7r°7r 7r°7r + channel. The left-hand side plots the 7r 7r 7r 7r + invariant 
mass distribution and the right-hand side is the 7r 7r°7r° invariant mass distribution. The 
continuous and dotted lines correspond to the old CLEO and new Novosibirsk current 
respectively. 




Figure 2: The z/ T 7r 7r°7r°7r + channel. The left-hand side plots the 7r°7r°7r + invariant mass 
distribution and the right-hand side is the tt°tt° invariant mass distribution. The continu- 
ous and dotted lines correspond to the old CLEO and new Novosibirsk current respectively. 
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7T°7T + 

Figure 3: The 7r 7T + invariant mass distribution for the u T 7i 7i°7i°7i + channel. The continu- 
ous and dotted lines correspond to the old CLEO and new Novosibirsk current respectively. 




Figure 4: The z/ r 7r + 7r + 7r~7r° channel. The left-hand side plots the 7r + 7r + 7r~7r° invariant 
mass distribution and the right-hand side is the 7r + 7r + 7r~ invariant mass distribution. The 
continuous and dotted lines correspond to the old CLEO and new Novosibirsk current 
respectively. 

implemented in the TAUOLA Monte Carlo library as options: 
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600 - 




7T + 7T + 7T 7Y + 7T 7v' 



Figure 5: The v T n + n + n~n channel. The left-hand side plots the 7r + 7r + 7r° invariant mass 
distribution and the right-hand side is the n + n^n invariant mass distribution. The con- 
tinuous and dotted lines correspond to the old CLEO and new Novosibirsk current respec- 
tively. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



7T 7T° 7T + TT + 

Figure 6: The v t 7t + 7t + 7t~tt channel. The left-hand side plots the 7r~ir° invariant mass dis- 
tribution and the right-hand side is the 7i + tt + invariant mass distribution. The continuous 
and dotted lines correspond to the old CLEO and new Novosibirsk current respectivevly. 
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0.8 
7T + 7T° 



Figure 7: The v T n + n + n~n G channel. The left-hand side plots the n + n~ invariant mass 
distribution and the right-hand side is the 7r + 7r° invariant mass distribution. The continu- 
ous and dotted lines correspond to the old CLEO and new Novosibirsk current respectively. 

• t + — > v t it + it it it channel 

In the first part of this chapter the invariant mass distributions for r + — > z/ T 7r + 7r 7r 7r° 
channel are shown. In Fig. 1 we show the invariant mass distribution for 7r 7r°7r°7r + (left- 
hand side plot) and 7r 7r 7r° (right-hand side plot) systems. In Fig. 2 we show the invariant 
mass distribution for 7r°7r°7r + (left-hand side plot) and 7r°7r° (right-hand side plot) systems, 
and in Fig. 3 the invariant mass distribution for 7r°7r + system. In all plots the continuous 
and dotted lines correspond to the old (1998) CLEO and new Novosibirsk current. 

• t + — > z/ r 7r + 7r _ 7r + 7r° channel 

In the second part of this chapter the invariant mass distributions for r + — > z/ r 7r + 7r _ 7r + 7r' 
channel are shown. In Fig. 4 we show the invariant mass distribution for 7r + 7r + 7r~7r° (left- 
hand side plot) and 7r + 7r + 7r~ (right-hand side plot) systems. In Fig. 5 we show the 
invariant mass distribution for 7r + 7r + 7r° (left-hand side plot) and 7r + 7r~7r° (right-hand side 
plot) systems. In Fig. 6 we show the invariant mass distribution for n + 7i° (left-hand side 
plot) and 7r + 7r + (right-hand side plot) systems. In Fig. 7 we show the invariant mass 
distribution for 7i + n^ (left-hand side plot) and 7r + 7r° (right-hand side plot) systems. 

In some cases agreement is only qualitatively correct, but one should have in mind 
rather limited data samples available at that time 4 . We expect the Novosibirsk model to 

4 The largest differences are present in the r+ — > 9 t tt + n~ tt + n a decay channel. They can be sub- 
stantially diminished, if the ujtt contribution to the current used in the CLEO model is appropriately 
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represent a substantial improvement over the old one. 

8 Summary 

The new parameterization of 4n form factors for the TAUOLA package is now available. 
The form factors are completely defined from the information in the paper. The partic- 
ular strength of the model used in their definition relies on its success in describing high 
statistics low energy e + e _ data at y/s < 1.4 GeV (future experiments at the upgraded 
collider VEPP-2000 will extend the energy range to the r lepton mass [27]). The CVC 
hypothesis has been instrumental in constructing predictions for the r decays. The cor- 
rectness of the code was checked by its comparison to the other generator based on the 
same data. Results obtained from our new form factors can be now used for comparisons 
with other models and r lepton decay data directly under conditions of any present or 
future experiment. 
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Appendix A 



Feynman diagrams for r + decays into 4tt 

Taking into account the quantum numbers of all intermediate and final states, see 
Table 3, the following Feynman diagrams can be written for the r + decay into z/ T 7r + 7r 7r 7r° 
and P T 7r + 7r~7r + 7r°, if ai(1260)7r and urn intermediate states are assumed. 

Table 3: The quantum numbers of mesons. 



Intermediate 






state 


I G 


jPC 


p(770) (p) 


1+ 


1— 




1" 


0" 


7T° 


1" 


o-+ 


ai(1260) 


1" 


1++ 


w(782) 


0" 


1— 


a 


0+ 


0++ 




Figure 8: Diagrams for the t + — > ^V^ + (9i)^ (Q , 2) 7r0 (Q , 3) 7 '" (9 , 4) decay via the p + — > ax7T 
p7T7T — > 7r7T7r7r intermediate states. 
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Figure 9: Diagrams for the r + — > ^T^iQi)^ ^)^ ^)^ ^) decay channel via the p 
a\K — > (J7T7T — > 7T7r7T7r intermediate states. 




Figure 10: Diagrams for the r + — > z/ r 7r + (g 1 )7r (f^) 71 "^^) 71 " !^) decay via the p + —* co>7r 
p7T7T — > 7r7T7T7r intermediate states. 
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Figure 12: Diagrams for the t + — > z/ r 7r + (gi)7r (<?2)7r + (<73)7r (<?4) deca?/ via the p 
(jTiTi — > 7T7r7r7r intermediate states 

Appendix B 



Tables of numerical values for functions G(Q 2 ) 
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